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formulate  a  comparable  capability  for  intermediate  conditions  such  as  grass 
covered  surfaces  where  both  the  terrain  surface  and  vegetation  influence  the 
average  surface  temperature. 


This  report  presents  a  procedure,  named  VEGIE,  that  predicts  the  tempera¬ 
ture  of  terrain  surfaces  which  contain  a  simple  layer  of  vegetation.  VEGIE  is 
designed  as  an  interim  procedure  for  immediate  application  in  lieu  of  more 
sophisticated  and  theoretical  treatments  of  this  problem.  Operational  flexi¬ 
bility  and  simplicity  are  preserved  by  using  VEGIE  as  a  submodel  to  the  Ter¬ 
rain  Surface  Temperature  Model  (TSTM)  developed  previously  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station.  The  TSTM  predicts  the  surface  tempera¬ 
ture  of  nonvegetated  layer  and  planar  features  using  material  thermophysical 

properties  and  meteorological  conditions.  ^ 

/ 

The  frameworks  of  the  TSTM  and  VEGIE  are  presented  along  with  the  strat¬ 
egy  for  their  joint  application  to  temperature  predictions  for  simply  vege¬ 
tated  surfaces.  The  sensitivity  of  the  temperature  predictions  to  the  inputs 
required  by  VEGIE  are  examined  to  assist  in  determining  the  limits  of  Its 
validity.  A  limited  validation  is  presented  using  measured  data  for  a  site 
in  West  Germany.  Results  from  VEGIE  are  also  compared  to  the  output  of  a  more 
sophisticated  (and  complete)  vegetation  temperature  model  developed  at  the 
Colorado  State  University  (CSU) . 

'>  These  models  are  applied  to  a  deciduous  and  a  coniferous  canopy  where 
the  CSU  model  is  valid  but  VEGIE  is  not.  VEGIE  is  applied  to  the  problem  of 
estimating  thermal  signatures  for  terrain  surfaces  with  less  than  total 
foliage  cover,  the  consequences  of  changes  in  foliage  cover,  differing  emis- 
sivltlea  for  the  soil  and  vegetation,  and  reflected  sky  radiation.. 

The  data  and  analyses  presented  herein  demonstrate  that  the  simplified 
treatment  of  energy  budgets  used  in  VEGIE  provide  realistic  results  for  the 
situations  that  VEGIE  was  designed  to  handle.  Under  some  conditions,  when 
sun  snd  view  geometry  are  important  or  whers  canopies  hsve  large  horizontal 
variations  to  outline  a  few,  the  results  from  VEGIE  may  be  of  limited  value. 
VEGIE  would  appear  to  perform  well  in  moderate  environmental  situations  such 
as  lawns,  pastures,  and  rangelands. 
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The  study  reported  herein  was  conducted  by  personnel  of  the  U.  S. 
Army  Engineer  Waterways  Experiment  Station  (WES)  from  1  October  1979  to 
1  October  1980.  The  study  was  done  under  Department  of  the  Army  Projects 
No.  4A762730AT42,  Task  A4,  Terrain/Operations  Simulation,  Work  Unit  003, 
Electromagnetic  Target  Surround  Characteristics  in  Natural  Terrains,  and 
No.  4A762719AT40,  Task  CO,  Theater  of  Operations  Construction,  Work 
Unit  006,  Fixed  Installation  Camouflage  Methods  and  Materials. 

The  study  was  conducted  under  the  general  supervision  of  Dr.  John 
Harrison,  Chief  of  the  Environmental  Laboratory,  and  Mr.  Bob  Benn,  Chief 
of  the  Environmental  Systems  Division,  and  under  the  direct  supervision 
of  Dr.  Lewis  E.  Link,  Jr.,  Chief  of  the  Environmental  Constraints  Group 
(ECG).  The  development  of  the  mathematical  model  presented  herein  was 
accomplished  primarily  by  Dr.  Lee  Balick,  on  assignment  to  ECG  from  Col¬ 
orado  State  University  (CSU).  Assistance  was  received  from  Messrs.  Randy 
Scoggins  and  Curt  Gladen,  ECG.  This  report  was  prepared  primarily  by 
Dr.  Balick  with  technical  assistance  from  Dr.  Link  and  Mr.  Scoggins, 
and  administrative  assistance  from  Ms.  Patti  Burke. 

Forest  canopy  thermal  signature  and  model  data  were  obtained  with 
the  collaboration  of  Dr.  Leo  J.  Fritschen  (College  of  Forest  Resources 
of  the  University  of  Washington,  Seattle,  Wash.),  Dr.  Boyd  A.  Hutchison 
(Atmospheric  Turbulence  and  Diffusion  Laboratory,  National  Oceanic  and 
Atmospheric  Administration,  Oak  Ridge,  Tenn.)  and  Dr,  James  A.  Smith 
(Department  of  Forest  and  Wood  Science,  CSU,  Fort  Collins,  Colo.).  In 
addition  to  their  scientific  collaboration,  Drs.  Fritschen  and  Hutchison 
were  hosts  of  field  experiments  at  their  research  sites  which  were  eoop- 
erative  efforts  between  their  organizations ,  CSU  and  WES .  Dr.  Smith 
led  the  development  of  the  CSU  thermal  model,  which  provided  temperature 
predictions  for  the  grass  and  forest  canopies. 

Commanders  and  Directors  of  the  WES  during  the  conduct  of  the  study 
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were  COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE.  Technical 
Director  was  Mr.  Fred  R.  Brown. 
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INCLUSION  OF  A  SIMPLE  VEGETATION  LAYER  IN  TERRAIN 


TEMPERATURE  MODELS  FOR  THERMAL  INFRARED 
(IR)  SIGNATURE  PREDICTION 

PART  I :  INTRODUCTION 

Background 

1.  Recent  analytical  and  theoretical  studies  of  the  thermal  in¬ 
frared  (IR)  emission  characteristics  of  terrain  surfaces  have  generally 
ignored  the  effects  of  vegetation  on  thermal  IR  signatures  (Gillespie 
and  Kahle  1977;  Pratt  and  Ellyett  1979;  Holmes,  Nuesch,  and  Vincent 
1980).  This  seems  to  be  due,  in  part,  to  a  lack  of  a  usable  tool  for 
examining  the  effects  of  vegetation  in  a  real  world.  In  truth  the  prob¬ 
lem  is  very  complex.  Models  of  vegetation  temperature  which  can  be 
directly  applied  to  remote  sensing  problems  (Kimes,  Smith,  and  Ranson 
1979;  Norman,  1979;  Smith  et  al.  1981a)  are  complex  and  require  careful 
specification  of  intracanopy  meteorological  conditions,  canopy  structure, 
and  biophysical  characteristics  which  are  not  often  available.  Observa¬ 
tional  and  analytical  studies  are  consistent  with  the  models  in  detailing 
the  complexity  of  the  problem  (Hiller  1971;  Heilman  et  al.  1976;  Bonn 
1978;  Byrne  et  al.  1979;  Kimes  1979;  Kimes  et  al.  1980;  Hillard  et  al. 
1980;  Soer  1980;  Balick  and  Wilson  1981). 

2.  Ongoing  work  at  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  has  as  a  major  goal  the  ability  to  realistically  predict 
thermal  signatures  of  natural  and  cultural  terrain  features  for  any 
time-of-day  and  weather  conditions.  Work  on  models  for  complex  vegeta¬ 
tion  canopies  and  planar  unvegetated  surfaces  has  pointed  to  the  need  to 
fill  the  gap  between  these  extremes.  A  model  for  predicting  the  time- 
varying  temperature  for  an  unvegttsted  planar  surface  using  material 
thermophysical  properties  and  meteorological  conditions  has  been  devel¬ 
oped  under  previous  WES  research  (Balick  et  al.  1981).  This  model,  the 
Terrain  Surface  Temperature  Hodel  (TSTH),  is  discussed  in  Part  II  of 
this  report.  Logically,  the  approach  taken  for  the  work  reported  herein, 
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the  next  step,  was  to  develop  a  module,  or  submodel,  for  use  in  con¬ 
junction  with  the  TSTM  to  account  for  the  dominant  effects  of  a  simple 
layer  of  vegetation  on  thermal  IR  signatures  of  the  terrain  surface. 
Such  a  capability  would  be  useful  even  if  it  only  applied  to  the  sim¬ 
plest  of  canopies  in  nonextreme  environments.  It  would  allow  the  TSTM 
to  be  extended  to  areas  of  lawn,  pasture,  and  perhaps  rangelands.  The 
module  developed  in  this  context  has  been  named  VEGIE  and  is  presented 
in  this  paper. 


Objectives  and  Scope 

3.  The  specific  objective  of  the  work  presented  herein  was  to  de¬ 
velop  the  capability  to  predict  the  temperature  of  terrain  scene  elements 
which  contain  a  simple  layer  of  vegetation  and  to  diagnose  the  effect  of 
vegetation  on  remotely  sensed  temperatures  of  terrain  elements.  More 
complete  and  theoretical  treatments  of  this  problem  are  under  way;  VEGIE 
is  designed  for  immediate  application  to  thermal  IR  signature  prediction 
and  analytical  studies.  It  was  also  required  that  the  operational  flex¬ 
ibility  and  simplicity  of  the  planar  surface  model  be  maintained.  Some 
of  the  more  important  of  these  characteristics  are: 

a.  Time  dependence  and  fast  response  to  environmental 
changes . 

b.  Air  temperature  considered  a  state  variable. 

c.  Materials  treated  as  horizontally  and  vertically 
homogeneous  layers. 

d.  Precipitation  and  condensation  not  considered. 

e.  Spectral  characteristics  not  considered  at  this  stage 
of  thermal  signature  prediction. 

f.  Sensible  and  latent  heat  transfers  included. 

g.  Cloud  type  and  amount  considered. 

Input  information  additional  to  that  used  in  the  TSTM  is  minimized  by 
the  use  of  empirical  and  quasi-eapirical  relationships  for  many  param¬ 
eters.  Most  of  these  relationships  use  coefficients  which  are  taken 
from  the  literature  and  are  not  unique.  Most  of  them  could  be  replaced 
by  measurements  if  measurements  are  available.  Therefore,  this  report 
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presents  details  of  only  the  more  basic  concepts  of  the  VEGIE  module. 

The  entire  set  of  equations  of  the  present  configuration  of  VEGIE  is, 
however,  given  in  Appendix  A  along  with  a  list  of  basic  symbols. 

4.  The  sensitivity  of  temperature  estimates  to  the  additional 
inputs  required  by  VEGIE  is  examined.  This  is  intended  to  examine  the 
behavior  of  VEGIE  and  to  help  determine  the  limits  of  its  validity.  (A 
qualitative  description  of  the  sensitivity  of  TSTM  is  given  along  with  a 
brief  description  of  the  model.)  The  model  is  then  applied  to  a  prob¬ 
lem  in  terrain  infrared  signature  prediction.  In  a  scene  element  with 
less  than  total  foliage  cover  there  exists  a  mixture  of  materials  (fo¬ 
liage,  soil)  with  emissivities  that  may  be  different.  The  consequences 
of  such  mixtures  are  examined  using  VEGIE,  with  and  without  considera¬ 
tion  of  reflected  sky  radiation. 

5.  Results  from  VEGIE  are  compared  to  two  day*  r>f  data  obtained 
in  Germany  in  order  to  obtain  a  partial  model  validation.  These  results 
are  also  compared  with  output  from  a  more  complete  vegetation  model  de¬ 
veloped  at  Colorado  State  University  (Smith  et  al.  1981a,  1981b).  Fi¬ 
nally,  interesting  results  from  applying  VEGIE  to  forests  are  presented. 
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PART  II:  MODEL  FRAMEWORK 


Overview 


6.  The  starting  point  for  the  development  of  the  vegetation  sub¬ 
model  was  a  temperature  prediction  scheme  developed  by  Deardorff  (1978). 
His  procedures  seem  to  be  designed  to  provide  an  efficient  way  to  account 
for  boundary  layer  transfer  of  sensible  and  latent  heat  in  atmospheric 
circulation  models.  Because  Deardorff 's  procedure  was  developed  for  a 
different  application  and  merged  with  different  computational  techniques, 
only  portions  of  his  work  were  found  to  be  directly  usable.  These  in¬ 
clude  the  energy  budget  equations  for  the  foliage  and  the  ground  surface, 
his  treatment  of  foliage  cover,  and  his  techniques  for  scaling  and  inter¬ 
polating  parameters  according  to  the  degree  of  vegetative  cover.  Also, 
his  equations  for  sensible  and  latent  heat  transfer  from  the  ground  were 
preserved.  Major  changes  were  aude  in  the  evaluation  of  terms  of  the 
foliage  energy  budget,  the  numerical  routine  for  solving  the  energy 
budgets,  and  the  overall  computational  sequence.  Still,  the  key  concepts 
are  derived  from  Deardorff 's  work. 

7.  The  flow  of  calculations  between  major  components  of  the  TSTM/ 
VEG1E  system  is  diagrammed  in  Figure  1.  At  block  l,  the  model  is  incre¬ 
mented  one  time  step  (or  initialized  or  terminated)  in  accordance  with 
procedures  established  for  the  TSTM.  The  flow  then  separates,  depending 
on  whether  or  not  VEGIE  is  to  be  used.  If  not,  a  nonvcgetated  surface 
energy  budget  is  evaluated  as  an  upper  boundary  condition  (block  3)  for 
the  solution  of  the  equation  for  b*«t  transfer  through  the  terrain  ma¬ 
terials  (block  6).  The  solution  for  the  surface  temperature  comes  from 
the  evaluation  of  the  surface  energy  budget  equation,  but  solving  the 
heat  transfer  equation  is  necessary  to  estimate  the  heat  conduction  term 
and  to  estimate  the  distribution  of  heat  in  the  terrain  materials.  If 
VEGIE  is  to  be  used,  the  energy  budget  of  block  3  is  replaced  by  blocks  4 
and  5.  Blocks  4  and  S  constitute  VEGIE.  Block  4  is  an  evaluation  of  the 
energy  budget  for  the  foliage  that  includes  a  contribution  from  the 
ground  surface.  Evaluation  of  the  ground  energy  budget  is  done  in 


7 


TERRAIN 

SURFACE 

temperature 

PREDICTION 


EVACUATE 
GROUND  SURFACE 
ENERGY  BUDGET  AND 
SOLVE  FOR  GROUND 
TEMPERATURE  IT#) 

< . I 


EVALUATE 

terrain  surface 

ENERGY  BUDGET 
AND  SOLVE  FOR  SUR¬ 
FACE  TEMPERATURE 


CALCULATE 
CONCLUSIVE  HEAT 
FLUXES  AND  TEMPS 
WTH1N  TERRAIN 
MATERIALS 


Fiture  1.  Sequence  of  calculation  cf  the  «a»jor 
components  of  the  TSTH/VEGIE  system 


block  5,  which  includes  a  contribution  from  the  foliage  layer.  (Solar 
and  thermal  IR  from  the  sky  incident  on  a  plane  above  the  vegetation 
are  determined  by  procedures  developed  for  the  TSTM  and  are  taken  here 
as  given.)  Some  details  of  the  calculation  of  sky  thermal  IR  radiation 
are  given  in  Appendix  A.  Solutions  of  temperature  for  the  foliage 
(from  block  4)  and  the  ground  (block  5)  are  performed  by  a  simple  root- 
findiug  algorithm  and  are  combined  according  to  the  proportion  of  fo¬ 
liage  cover  to  yield  an  average  or  effective  temperature  of  the  vegetated 
surface.  The  ground  energy  budget  is  then  used  in  the  evaluation  of  heat 
flow  in  the  terrain  (block  6),  and  the  program  returns  to  block  1.  The 
calculations  of  blocks  4  and  5  are  the  primary  subjects  of  this  report 
but  a  brief  description  of  the  TSTM  is  warranted. 

TSTM 


8.  The  TSTM  solves  the  one-dimensional  time-dependent  heat  flow 
equation  for  material  systems  of  one  to  six  horizontally  and  vertically 
homogeneous  layers.  The  equation  is  written  as 

Sllt.t)  „  ,(l)  ii&n  „) 

9t  az2 

where  T  is  temperature,  z  in  distance  into  the  material  system, 
a  is  the  thermal  diffusivity,  and  t  is  time.  The  solution  of  this 
equation  is  done  with  straight  explicit  numerical  techniques  given  an 
initial  temperature  profile  and  upper  and  lower  boundary  conditions 
which  are  energy  budgets.  The  lower  boundary  condition  is  user- 
specified  in  a  very  simple  fashion  with  three  available  options.  The 
upper  boundary  is  of  greater  interest;  from  it  is  derived  the  surface 
temperature  estimate. 

9.  The  heat  budget  of  the  top  surface  can  be  written  as  the  sum 
of  energy  flux  density  components: 

S  +  R  ,  +  R  +  H  +  E  +  G  =  0  (2) 

8*  T 
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where 

S  is  the  solar  radiation  absorbed  at  the  surface 

R  ,  is  the  downwelling  thermal  IR  energy  from  the  sky 
8* 

R^  is  the  thermal  IR  energy  emitted  by  the  surface 
H  is  the  sensible  heat  exchange 
E  is  the  evaporative  heat  exchange 

G  is  the  conductive  transfer  between  the  surface  and  the  terrain 
material 

This  equation  is  rewritten  as  a  nonlinear  function,  F  ,  of  the  surface 

temperature  such  that  the  updated  surface  temperature  is  a  zero  of  F  . 

The  Newton-Raphson  algorithm  is  used  to  locate  the  value  of  T  such 

s 

that  F  is  zero. 

10.  Each  term  of  Equation  2  is  evaluated  separately.  S  is  the 
global  solar  insolation  multiplied  by  the  absorptivity.  Insolation  is 
usually  an  input,  but  it  is  estimated  by  TSTM  if  it  is  not  in  the  input 
file.  Downwelling  thermal  IR  from  the  sky  is  estimated  with  the  Brunt 
equation  with  adjustments  for  cloud  type  and  amount  (Sellers  1965  or 
Oke  1978).  Rf  is  the  graybody  radiant  emittance  and  H  and  E  are 
evaluated  with  the  "aerodynamic"  approach  to  turbulent  heat  transfer 
(Oke  1978)  with  adjustments  for  stability.  G  is  approximated  by 


G  =  K 


(3) 


where  K  is  the  heat  conductivity  of  the  top  layer,  Ts_j  is  the  temper¬ 
ature  of  the  first  grid  point  below  the  surface,  and  Az  is  the  distance 
between  that  point  and  the  surface. 

11.  Inputs  to  TSTM  include  atmospheric  constants  and  time-varying 
atmospheric  and  material  properties.  Atmospheric  constants  required  are 
atmospheric  pressure  and  instrument  shelter  height  above  the  surface. 

The  shelter  height  value  represents  the  height  above  the  ground  that  air 
temperature  and  wind  speed  are  measured.  Atmospheric  time-varying  data 
required  include  air  temperature,  relative  humidity,  cloud  cover,  wind 
speed,  and  total  insolation.  Solar  insolation  can  also  be  confuted  as 
previously  mentioned.  Material  properties  are  needed  for  the  surface 
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and  each  of  the  layers.  Surface  properties  required  are  thermal  erais- 
sivity,  optical  absorptivity,  and  percent  saturation  of  the  surface. 

Each  layer  is  defined  by  its  thickness,  thermal  diffusivity,  and  heat 
conductivity. 

12.  Table  1  presents  a  qualitative  assessment  of  model  sensitivity 
to  its  input  data.  Figure  2  illustrates  the  response  of  the  model  to 
changes  of  time-varying  atmospheric  data;  in  this  case  the  data  are 
hourly  and  the  only  variable  changes  are  cloud  cover.  TSTM  is  described 
in  detail  by  Balick  et  al.  (1983).  A  conceptually  similar  model  has 
been  developed  by  Kakle  (1977). 

13.  Operationally,  the  only  change  needed  to  run  VEGIE  is  to  add 
a  single  line  of  data  to  the  input  file.  If  these  data  are  present,  the 
program  bypasses  the  above  top  boundary  energy  calculations  and  proceeds 
as  described  below. 

VEGIE 

14.  Description  of  the  VEGIE  module  is  done  in  four  parts;  the 
foliage  energy  budget  (Figure  1,  block  4),  the  ground  energy  budget 
(Figure  1,  block  5),  the  root-finding  algorithm  used  to  solve  each  energy 
budget  for  temperature,  and  the  combination  of  these  temperatures  to  form 
an  average  temperature  (Figure  1,  block  7).  The  overall  sequence  of  cal¬ 
culations  has  been  discussed.  The  complete  set  of  equations  is  given  in 
Appendix  A  and  their  sources  are  referenced  there. 

Foliage  energy  budget 

15.  The  equation  for  the  foliage  energy  budget,  ,  is  adopted 
directly  from  Deardorff  (1978)  and  is 

Ff  =  OfC«fS  ♦  +  Rfl)  -  Hf  -  Ef  (4) 

where  o f  is  the  fraction  of  foliage  cover,  a  is  the  shortwave  (solar) 
absorptivity,  c  is  graybody  emissivity,  the  subscript  f  pertains  to 
the  foliage,  and  Rft  is  a  combined  "net"  thermal  IR  term  for  interaction 
between  the  foliage,  ground,  and  their  loss  to  the  sky.  The  equation 
for  Rq  is  in  Appendix  A  and  is  a  function  of  the  foliage  and  ground 
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Table  1 


Relative  Response  of  TSTM  to  Variation  of  Input  Parameters 


Very  Sensitive 


Air  temperature 

Solar  absorption 

Thermal  emissivity 

Initial  temperature 
profile 

Saturation 

Cloud  cover  (low 
level  clouds) 

Top  layer  heat 
conductivity 


Moderately  Sensitive 


Relative  humidity 
Shelter  height 


Wind  speed 

Cloud  cover  (middle 
lenel  clouds) 

Cloud  type: 

(between  high,  middle 
and  low  level  clouds) 


Insensitive 


Air  pressure 

Cloud  cover  (high 
level  clouds) 

Thermal  diffusivity 

Time  step* 

Grid  spacing* 

24-hr  repetitions 

Bottom  boundary 
flux** 


*  Mot  sensitive,  provided  model  is  numerically  stable. 

**  Not  sensitive  for  thick  systems  with  relatively  low  beat 
conductivity  Including  most  soils. 
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CASE  1  NO  CLOUDS  ALL  HOURS 
CASE  2  NO  CLOUDS  EXCEPT  OVERCAST  1100-1300  HR 
CASE  3  NO  CLOUDS  EXCEPT  OVERCAST  0800-1300  HR 
CASE  4  NO  CLOUDS  EXCEPT  OVERCAST  0100-0600  HR 


_ I _ _ _ I _ I _ I  i _ I _ L 

0  4  8  12  16  20  24 

TIME  OF  OAY,  HOURS 

Figure  2.  TSTM  predictions  of  road  surface  temperature 
with  time-varying  cloud  cover  conditions  (from  Balick 

et  al.  1981) 


emissivities  (e^  ,  )  and  temperatures  (T^  ,  T  ).  The  parameters 

af  ,  a^,  ,  ef  ,  are  inputs  to  the  model,  S  and  Rg^  are  inputs 

to  and/or  calculated  by  TSTM,  T  is  the  ground  temperature  estimate 
from  the  previous  time  step,  and  is  determined  by  the  root-finding 

algorithm.  Sensible  heat  transfers  from  the  foliage  are  primarily  func¬ 
tions  of  temperature,  wind  speed,  and  vapor  pressure  differences  between 
the  foliage  or  the  air  adjacent  to  it  and  the  air  at  the  instrument 
shelter  height.  Energy  storage  and  conduction  by  the  foliage  are 
neglected. 

16.  Values  of  five  foliage  parameters  are  required  for  VEGIE.  They 


are: 


a. 

V 

foliage  cover  fraction,  0  <  <  1 

b. 

X: 

state  of  vegetation,  x  >  0 

c. 

V 

graybody  emissivity,  0  <  ef  <  1 

d. 

V 

shortwave  (solar)  absorptivity  of  the  layer, 

0  <  af  <  1 

e. 

V 

foliage  height,  >  0 

17.  Deardorff  defines  crf  as  an  area  average  shielding  factor 
associated  with  the  degree  to  which  foliage  prevents  shortwave  radiation 
from  reaching  the  ground.  It  is  not  the  same  as  a  visual  cover  because 
of  single  or  multiple  shortwave  scattering  by  the  leaves;  is  less 
than  the  degree  visual  cover  to  some  extent.  The  range  of  a{  is  0  to 
1  and  is  assumed  not  to  vary  with  sun  zenith  or  azimuth  angle.  Deardorff 
provides  typical  values  from  Geiger  (1965)  as  follows:  0.82  for  meadow 
grass,  0.95  for  30-cm-high  clover,  0.83  for  80-cm-high  winter  rye,  0.30 
for  summer  barley  12-15-cm  high,  and  from  0.4  to  0.98  for  various  stands 
of  trees.  The  foliage  cover  parameter  can  be  roughly  related  to  leaf 
area  index  (LAI)  for  ground  vegetation  by: 


Of  4  LAI/7  (5) 

X  is  used  as  a  multiplier  of  the  stomatal  resistance  function  in  VEGIE. 
An  empirical  relationship  of  stomatal  resistance  as  a  function  of  inso¬ 
lation  for  typical  .unstressed  plants  is  used  in  VEGIE.  Use  of  %  then 
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allows  rough  adjustment  for  deviations  from  this  curve  by  the  user. 
Normally  x  =  1.0  is  used,  but  other  values  can  be  chosen  to  adjust 
stomatal  resistance  for  moisture  stress,  senescence,  or  other  factors. 
The  meaning  of  the  other  parameters  is  clear. 

18.  All  of  the  radiant  energy  terms  of  the  foliage  energy  budget, 
Equation  4,  are  weighted  by  ,  but  the  turbulent  transfer  terms  are 
not  weighted.  These  terms  are  evaluated  for  a  unit  of  ground  area. 

(One  slight  exception  is  in  determining  the  value  for  the  canopy 
resistance  to  water  vapor  diffusion  for  the  term;  see  Appendix  A.) 

Equation  4  is  solved  for  the  value  of  T^  that  makes  =  0  . 

Ground  energy  budget 

19.  Following  Deardorff  again,  the  energy  budget  for  the  ground 
is  written  as 


F  =  (1  -  oJa  S  +  R-R-H  -E  -  G  (6) 

8  f  8  g*  8+  8  8 

The  adjustment  of  insolation  for  foliage  cover  is  clear  and  no  adjustment 

is  appropriate  for  the  conduction  term,  G  .  Otherwise,  adjustments  for 

foliage  density  are  made  in  evaluating  the  terms.  Terms  for  incoming 

and  outgoing  thermal  IR  flux  density,  R^  and  R^t  ,  are  each  functions 

of  T  ,  T.  ,  e  ,  ef  ,  and  R  .  and  equations  include  reflection, 
gig*  s* 

Evaluation  of  the  turbulent  energy  transfer  terms  for  the  foliage  energy 
budget  is  done  per  unit  ground  area  and  with  the  assumption  that  the  ex¬ 
change  takes  place  directly  between  the  foliage  and  the  air  above  it. 
Turbulent  transfer  at  the  ground  is  still  per  unit  area  but  is  buffered 
by  the  layer  of  vegetation  and,  thus,  is  a  function  of  a ^  .  Addition¬ 
ally,  the  transfer  occurs  between  the  ground  and  the  air  within  the  fo¬ 
liage  layer  which  is  an  arbitrary  mixture  of  conditions  at  the  ground, 
foliage,  and  air  at  shelter  height.  The  evaluation  of  these  terms  fol¬ 
lows  Deardorff  very  closely.  Conduction  is  estimated  with  Equation  3  as 
in  the  base  ground  energy  budget.  The  root-finding  algorithm  is  used  to 

find  the  value  of  T  that  makes  F  =  0  . 

•  8 

Root-finding  algorithm 

20.  The  objective  of  the  root-finding  algorithm  is  to  find  the 
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temperature,  T^  or  T  ,  for  which  the  sum  of  the  energy  budget  com¬ 
ponents  is  zero  (or  within  a  specified  range  from  zero).  Solved  for 
temperature,  the  energy  budgets  are  fourth-order  equations  but  experience 
has  shown  that  there  is  only  one  real  and  reasonable  root  given  realistic 
conditions.  Therefore,  a  very  simple  and  efficient  algorithm,  similar 
to  those  in  software  packages  used  for  many  microcomputers,  was  adapted 
for  VEGIE  and  is  called  the  regula  falsi  technique  (Scheid  1969). 

21.  Figure  3  schematically  illustrates  the  steps  in  the  algorithm. 
An  initial  guess  of  temperature,  T  ,  is  made  which  is  unrealistically 
low  (say  200°K)  such  as  point  A,  in  Figure  3a.  The  energy  budget, 

F  =  f(T)  ,  (Equation  6  or  8)  is  evaluated  there  and  its  sign  is  deter¬ 
mined.  F  is  reevaluated  at  progressively  higher  temperatures  at  regular 
steps  (say  5°K)  until  the  sign  of  F  changes.  Points  A  and  B  in  Fig¬ 
ure  3b  are  determined.  The  intercept  of  a  straight  line  between  points 
f(A)  and  f(B)  is  found;  point  C  in  Figure  3c.  Then  f(A)  becomes  f(C)  , 
a  new  line  between  f(A)  and  f(B)  is  determined,  a  new  intercept  is 
found  (Figure  3d)  and  so  on  until  F  is  less  than  some  assigned  value 
(F  <  0.001  for  VEGIE).  This  algorithm  is  simple,  straightforward,  and 
easily  adaptable  to  microcomputers  for  functions  that  behave  like  F  . 
Effective  temperature 

22.  The  simplest  procedure  for  combining  ground  and  vegetation 
temperature  would  be  a  simple  average  weighted  by  the  foliage  cover: 

f  =  ofTf  +  (1  -  of)Tg  (7) 

However,  when  temperatures  are  observed  radiometrically,  it  is  more 
appropriate  to  mix  the  radiant  exitance  from  the  two  materials  and 
solve  for  temperature  as  follows: 

[>  ,*10.25 

°fefTf  *  (1  *  °f)cgTgJ  (8) 

Except  as  otherwise  stated,  all  model  results  presented  in  following  sec¬ 
tions  are  effective  temperatures  as  defined  in  Equation  8.  The  effective 
temperature  value  is  the  primary  product  of  the  TSTN/VEGIE  system. 
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PART  III:  SENSITIVITY  AND  VERIFICATION 


23.  This  section  presents  a  brief  analysis  of  the  sensitivity  of 
temperature  estimates  to  the  input  vegetation  parameters  required  by 
VEGIE.  Model  estimates  are  then  compared  to  two  days  of  real  data  and 
estimates  from  a  more  complex  model  for  the  same  two  days.  Model  esti¬ 
mates  for  two  tree  canopies  are  presented  as  cases  where  the  concepts  of 
VEGIE  are  not  valid,  but  the  counterexamples  are  informative.  Lastly, 
VEGIE  is  applied  to  examining  the  effects  of  emissivity  mixtures  within 
scene  elements  on  thermal  IR  signature  prediction  and  analysis. 

Sensitivity 

24.  The  sensitivity  of  VEGIE  to  its  input  variables  was  examined 
by  allowing  each  of  the  variables  to  vary  over  a  large  range,  while  hold¬ 
ing  other  variables  the  same  for  one  arbitrary  but  real  diurnal  cycle  of 
environmental  conditions.  The  change  of  model  estimates  as  a  function  of 
parameter  value  is  then  isolated  for  that  case.  Figure  4  shows  the  sensi¬ 
tivity  of  model  estimates  to  changes  of  .  The  curves  are  the  highest 
and  lowest  temperature  for  a  complete  diurnal  cycle  at  different  foliage 
covers.  In  Figure  4a,  minimum  effective  temperature  estimates  are  almost 
unchanged  across  the  range  of  a ^  .  Maximum  effective  temperatures  are 
strongly  affected  by  changes  of  at  low  values  but  when  >  0.7 
there  is  almost  no  effect.  These  curves  are  a  combination  of  changes  in 
energy  transfers  and  the  mixture  of  ground  and  foliage  used  in  averaging 
to  obtain  effective  temperatures.  Similar  curves,  but  for  ground  and 
foliage  temperatures  separately,  are  given  in  Figure  4b.  This  removes 
the  averaging  ambiguity  and  clearly  illustrates  that  the  ground  tempera¬ 
ture  estimates  are  far  more  sensitive  to  changes  of  than  foliage 
temperature  estimates.  This  implies  that  foliage  temperatures  do  not 
change  auich  with  the  density  of  foliage  cover.  Expectedly,  maximum 
ground  temperature  estimates  are  more  affected  than  minimum  temperatures. 

25.  VEGIE  was  found  to  be  highly  insensitive  to  changes  of  the 
other  four  foliage  parameters.  For  one  parameter  at  a  time,  each  of  the 
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four  inputs  were  varied  over  an  unrealistically  large  range  of  values, 
and  neither  maximum  nor  minimum  temperature  estimates  changed  by  more 
than  1°C.  A  close  look  at  the  energy  budget  components  showed  that 
there  is  an  underplay  or  feedback  between  them.  That  is,  when  x  was 
raised  to  1000,  was  reduced  to  near  zero.  However,  other  terms, 

most  notably  ,  changed  to  compensate  for  the  reduction  of  latent 
heat  loss.  Such  feedbacks  do  exist  in  nature,  and  in  moderate  and  un¬ 
stressed  conditions,  temperatures  are  closely  coupled  to  air  temperature. 
The  results  with  VEGIE  are  consistent  with  this  situation.  Only  exten¬ 
sive  comparison  with  careful  measurements  can  indicate  the  extent  to 
which  these  model  responses  represent  real  conditions. 

Verification 

26.  Figure  S  contains  the  model  estimates  of  terrain  temperature 
and  air  temperature  over  two  diurnal  cycles;  1  September  1979  in  Fig¬ 
ure  5a  and  10  October  1979  in  Figure  5b.  Ranges  of  radiometric  tempera¬ 
tures  are  shown  as  vertical  bars  at  the  time  of  observation.  The  data 
were  obtained  at  Zweibruecken  Air  Force  Base,  West  Germany,  during  the 
course  of  a  North  Atlantic  Treaty  Organization  (NATO)  field  experiment. 

On  1  September  1979  the  weather  was  warm  with  the  sky  changing  from 
clear  to  0.4  or  0.5  cover  of  cirrus  clouds  after  OiJOO  local  standard 
time  (LST).  On  11  October  1979  the  weather  was  cooler  with  the  sky 
changing  from  clear  to  0.9  or  1.0  cover  of  stratocumulus  clouds  after 
about  0400  LST.  TSTM  estimates  of  solar  insolation  are  used  in  lieu  of 
observations  for  both  days.  Vegetation  is  rough,  unstressed  grass  about 
10  cm  high. 

27.  Examination  of  Figure  S  reveals  that  model  estimates  follow 
air  temperatures  closely  but  that  they  do  deviate  in  a  manner  that  would 
be  expected  giver  time  of  day  and  cloud  cover.  Hore  importantly,  these 
deviations  are  generally  of  the  same  signs  and  magnitudes  as  measure¬ 
ments  indicate  are  real.  The  measurements  on  1  September  1979  were  made 
by  a  variety  of  observers  using  different  instruments  at  different  sites. 
Thus,  observations  that  deviate  markedly  from  air  temperature  might  be 
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b.  ZWElBAOCKEN  AFB.  GE.  11  OCT  79 

Figure  S.  Model  estimates  end  radiometric  observations  o f  terrain 
temperatures  and  air  temperatures  for  two  days  at  Zweibrueckeu  Air 
Force  Base,  Vest  Germany,  a{  -  0.85  ,  x  -  1.0  ,  tf  =  0.98  , 

-  0.75  ,  »  10  cm 
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invalid.  If  so,  the  performance  of  VEGIE  must  be  considered  quite  good 
for  the  very  moderate  conditions  on  these  test  days. 

28.  What  accuracy  is  lost  in  using  a  simple  model  instead  of  a 
complex  model  on  the  grass  canopies  for  these  days?  Are  the  simplifica¬ 
tions  of  more  elaborate  theoretical  procedures  valid?  A  means  to  answer 
these  questions  more  completely  is  the  Thermal  Vegetation  Canopy  Model 
developed  at  Colorado  State  University  (CSU)  (Smith  et  al.  1981a,  b). 

This  model  was  run  on  the  same  canopy  and  for  the  same  two  days  as  in 
Figure  5.  Additionally,  the  same  TSTM  insolation  estimates  and  VEGIE 
ground  surface  temperature  estimates  were  used  in  the  CSU  model  calcula¬ 
tions.  (Ground  temperatures  are  a  required  input  to  the  CSU  model.) 

This  last  feature  implies  that  the  differences  in  results  are  due  to  dif¬ 
ferences  in  foliage  temperature  estimates  and  procedures  used  to  derive 
an  effective  temperature. 

29.  Estimates  from  the  two  models  are  very  similar  as  can  be  seen 
in  Figure  6.  Only  occasionally  are  the  results  from  the  two  models  more 
than  1°C  different,  mostly  for  the  morning  of  1  September  1979.  VEGIE 
estimates  of  the  maximum  temperature  on  both  days  are  about  1°C  warmer 
than  the  CSU  model  but  the  difference  is  short-lived.  The  source  of 
different  results  of  the  two  models  has  not  been  specifically  deter¬ 
mined.  In  any  case,  not  much  predictive  accuracy  is  lost  in  using  VEGIE 
on  the  simple  layer  of  vegetation;  much  time  and  effort  can  be  saved. 

30.  During  the  sussser  of  1979.  two  collaborative  experiments  were 
done  to  obtain  data  to  derive  and  validate  canopy  thermal  IR  signature 
models  in  forest  situations  (Smith  et  al.  1981a).  Two  complete  diurnal 
cycles  of  data  were  obtained  in  a  Douglas-fir  canopy  (28  m  high,  Cedar 
River,  near  Issaquah,  Wash.)  and  two  more  in  a  mixed  deciduous  canopy 
(21. S  m  high,  Walker  Branch  near  Oak  Ridge,  Tenn.).  The  temptation  to 
run  VEGIE  on  these  data  sets  was  not  resisted.  Forest  esnopies  are,  how¬ 
ever,  outside  of  the  range  of  conditions  intended  for  the  use  of  VEGIE. 

31.  Figure  7  contains  the  results  from  the  Wslker  Branch,  Tenn., 
deciduous  canopy  on  18  and  19  August  1979.  Unlike  the  estimates  for 
Germany,  ground  temperatures  for  the  CSU  model  were  measured  with  s  hand¬ 
held  radiometer.  The  two  models  agreed  very  well  at  night,  hut  VEGIE 
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ZWEIBRUCKEN  AFB,  GE,  11  OCT  79 

Figure  6.  Comparison  of  VEGIF.  and  Colorado  State  University 
model  estimates  and  observations  (VEGIE  estimates  and 
observations  are  the  same  as  in  Figure  5) 
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MIXED  DECIDUOUS  CANOPY 


Figure  7.  Comparison  of  models  and  observations  when  applied  to  a  21.5-m-tall 
mixed  deciduous  forest  canopy  at  Walker  Branch,  Tenn.,  on  18  and  19  August  1979 


seriously  overpredicts  temperature  during  the  day.  Similar  results  were 
obtained  for  the  Douglas-fir  canopy  at  Cedar  River,  Wash.,  on  4  and 
5  August  1979.  These  results  are  shown  in  Figure  8.  VEGIE  seems  to 
track  well  with  the  radiometric  measurements  made  from  the  ground,  but 
this  is  probably  not  significant.  Thermocouple  measurements  are  supe¬ 
rior  in  this  case  and  should  be  used  as  the  standard  for  comparison. 

Also  notable  is  a  tendency  for  VEGIE  to  predict  slightly  lower  temper¬ 
atures  than  the  CSU  model. 

32.  Close  examination  of  the  daytime  overestimation  of  forest  can¬ 
opy  temperatures  indicates  that  the  magnitude  of  the  error  is  more 
closely  related  to  insolation  than  any  other  environmental  variable. 

The  equations  in  VEGIE  describe  the  solar  irradiation  on  the  foliage  as 
occurring  on  a  porous  but  single  flat  surface.  In  nature  and  in  the  CSU 
model  (and  many  other  models  as  well),  the  solar  energy  is  vertically 
distributed  throughout  the  canopy.  That  simplification  for  VEGIE  may, 
at  least  in  part,  explain  VEGIE' s  tendency  to  overestimate  temperatures 
during  the  day.  The  lack  of  a  vertically  distributed  canopy  in  VEGIE 
may  also  be  responsible  for  its  nocturnal  underestimates.  It  would  be 
fortunate  if  these  speculations  were  accurate  because  the  problem  would 
be  greatly  reduced  for  grassy  canopies;  heights  are  on  the  order  of  20  cm 
instead  of  20  m.  There  is  only  a  hint  of  such  behavior  in  the  grass 
canopy  simulations  given  in  Figure  5. 

Application:  Mixed  Emissivities 

33.  Transformation  of  temperature  estimates  to  thermal  IR  signa¬ 
tures  requires  specification  of  euissivity.  When  a  simple  foliage  layer 
partially  covers  a  soil  surface,  a  sensor  receives  energy  emitted  from 
both  in  proportion  to  the  amount  of  cover,  their  emissivities,  and  their 
temperatures.  Additionally,  the  sensor  receives  thermal  IR  energy  re¬ 
flected  by  the  surface  materials  from  the  sky  and  surroundings.  Radia¬ 
tion  from  the  sky  and  this  reflection  varies  mainly  with  atmospheric 
temperature  and  water  vapor  and  clouds  and  is  then  time-dependent.  Vary¬ 
ing  mixtures  of  emissivity  also  affect  energy  flumes  and  transformations 
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to  some  extent.  Clearly,  the  effects  of  emissivity  mixtures  within  a 
scene  element  can  be  a  critical  issue  on  signature  prediction  and  analy¬ 
sis.  They  are  examined  in  two  ways  in  this  section.  First,  a  single 
arbitrary  case  is  presented  where  only  the  foliage  cover  and  ground 
emissivity  are  varied.  Secondly,  VEGIE  is  used  to  examine  the  more  com¬ 
plex  case  where  environmental  conditions  and  energy  budget  changes  are 
considered. 

34.  The  thermal  IR  energy  emitted  per  unit  time  and  area,  , 
for  a  surface  composed  of  foliage  and  ground  is 


Wt  =  af£foT*  +  (1 


CT-)£  OT^ 
f  g  g 


(9) 


For  a  graybody  near  thermal  equilibrium,  thermal  IR  reflectivity  is 
(1  -  e)  .  Thus,  a  sensor  pointed  down  at  a  grass/ground  surface  would 
receive  emitted  and  reflected  radiation  from  the  surface.  Ignoring 
multiple  reflection  and  atmospheric,  spectral,  and  directional  effects, 
this  quantity  is 


Wt  *  af6foTf  +  (1  ~  V^g 


+  (1  -  °f)(l  »  eg)eaoT*  +  °fU  ’  ef)eaoT J  (10) 


(For  a  more  complete  treatment  see  Lorenz  (1966))  where  the  subscript  a 
indicates  values  for  the  atmosphere.  from  Equations  9  and  10  can 
be  converted  to  average  effective  temperature  with 

/W.\0,25 

i=(;r)  OD 


Foliage  emissivities  are  generally  greater  than  0.95  and  usually  about 
0.98  (Cates  1980);  whereas,  mineral  soil  emissivities  can  be  considerably 
lower  (Taylor  1979;  Holmes,  Nueach,  and  Vincent  1980;  Buettner  and  Kern 
1966)  and,  in  the  field,  are  less  well  known.  Figure  9a  shows  the  varia¬ 
tion  of  T  with  changes  of  and  eg  for  an  arbitrary  case,  with  and 
without  consideration  of  reflection.  Figure  9b  shows  the  error  from 
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REFLECTtON  ERROR.  °C  TEMPERATURE. 


a. 


Figure  9.  Comparison  of  effective  temperatures  as  a  function 
of  foliage  cover  with  different  soil  emissivities  and  with 
and  without  considerations  for  reflected  sky  thermal  IR 

irradiance 
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neglecting  reflection  as  a  function  of  a ^  for  two  values  of  £ . 

When  reflection  is  not  considered,  an  uncertainty  of  0.04  in  £^  results 

in  a  4°C  error  in  T  at  =  0  and  about  1.5°C  at  0^  =  0.5  and  is 

small  above  =  0.8  .  Errors  of  emissivity  affect  the  reflection  terms 

in  the  opposite  direction  which  greatly  reduces  the  net  effect.  Errors 

caused  by  ignoring  reflection  in  this  case  are  greater  than  5°C  for 

e  <0.92  at  a  =  0  ,  greater  than  3°C  at  of  =  0.5  ,  and  always 

greater  than  1°C.  Treating  bare  ground  as  a  blackbody  (£  =  1  and  no 

8 

reflection)  results  in  an  error  of  8.5°C.  These  figures  are  for  a  rea¬ 
sonable  but  arbitrary  case. 

35.  The  use  of  VEGIE  easily  permits  a  similar  analysis  but  for 

the  course  of  an  entire  day  and  for  specific  environmental  and  terrain 
situations.  Since  reflection  is  a  function  of  atmospheric  temperature 
and  humidity  through  T  and  £  ,  it  is  subject  to  changes  in  time  as 

well  as  surface  composition.  The  use  of  VEGIE  to  diagnose  these  effects 
can  give  a  somewhat  different  picture  than  that  in  Figure  9  and  can  be 
very  important  in  efforts  to  use  and  predict  thermal  IR  signatures. 

36.  The  curves  plotted  in  Figure  10  are  model  estimates  of  maximum 
and  minimum  temperatures  as  a  function  of  0f  for  a  real  case.  The  case 
is  the  same  as  in  Figure  5a,  1  September  1979  at  Zweibruecken  Air  Force 
Base,  West  Germany.  Temperatures  are  defined  by  Equations  9,  10,  and  11. 
The  magnitude  of  the  reflection  effects  is  similar  to  those  in  Figure  9a 
but  is  larger  for  the  maximum  temperatures  than  the  minimum.  The  pre¬ 
dicted  range  of  T  as  a  function  of  0^  varies  from  about  1°C  for  min¬ 
imum  temperatures  to  about  9°C  for  maximum  temperatures;  both  extremes 
occur  with  reflection  considered.  Importantly,  the  shapes  of  the  curves 
in  Figures  9a  and  10  are  radically  different.  According  to  VEGIE,  the 
minimum  diurnal  range  is  found  at  intermediate  values  of  foliage  cover-- 
at  *  0.6  without  reflection  and  at  =  0.8  with  reflection. 

These  kinds  and  magnitudes  of  effects,  due  to  reflection  and  uncertainty 
of  effective  emissivity,  are  significant  for  the  prediction,  interpreta¬ 
tion,  and  analysis  of  thermal  IR  measurements. 
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and  ground  emiaaivity  (eg)  with  and  without  reflection  of 
sky  thermal  IR  irradiance 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


37.  The  highly  simplified  treatment  of  the  energy  budgets  of  fo¬ 
liage  and  ground  in  VEGIE  seems  valid  for  the  simple  cases  it  is  designed 
to  handle.  This  conclusion  is  indicated  by  the  results  obtained  for  the 
Zweibruecken  Air  Force  Base  grass  canopy,  as  discussed  in  paragraphs  26 
and  27.  However,  the  module  is  not  designed  to  be  valid  for  all  canopies. 
Results  from  the  application  of  VEGIE  to  forest  canopies,  presented  in 
paragraphs  30  and  31,  show  that  VEGIE  does  not  perform  well  in  these  com¬ 
plex  canopies.  VEGIE  estimates  of  temperature  are  sensitive  to  changes 
of  foliage  cover  but  not  very  sensitive  to  the  other  canopy  descriptors 
(paragraphs  24  and  25).  This  implies  that  VEGIE  may  not  perform  well 
under  environmental  conditions  where  turbulent  heat  transfers  and  plant 
responses  are  inadequately  described. 

38.  It  is  recommended  that  additional  sets  of  validation  or  test 
data  be  obtained  in  order  to  better  define  the  canopy  and  environmental 
conditions  where  VEGIE  is  valid.  Additional  but  limited  development  of 
simple  models  seems  warranted.  It  is  recommended  that  a  simple  treat¬ 
ment  of  canopy  radiation  penetration  be  developed  for  VEGIE  to  account 
for  the  vertical  distribution  of  foliage.  The  long-term  development  of 
simple  models,  however,  is  limited  for  several  reasons.  These  include 
the  necessity  to  adequately  specify  intracanopy  environmental  conditions, 
plant  physiological  parameters,  and  the  structure  of  the  canopy  itself. 
More  importantly,  there  is  a  lack  of  fundamental  understanding  of  turbu¬ 
lent  transfers  of  energy  in  plant  canopies;  better  predictive  capability 
awaits  a  better  theoretical  base.  Long-range  advances  in  operational 
prediction  of  thermal  IR  signature  of  vegetation  will  likely  best  be 
made  through  the  simplification  and  adaptation  of  complex  models  and 
theoretical  developments. 
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APPENDIX  A:  ENERGY  BUDGETS  AND  TERM  DEFINITIONS 


Foliage 


Ff  ■  +  V  -  "f  •  Ef 

cr ^.S :  a{  is  input;  S  is  input  or  estimated  in  TSTM 

£fRg+:  ef  is  input;  Rg+  from  TSTM  (Balick  et  al.  1981*) 

R  =  R.  -  R0  where 
n  1  2 

Ri *  (efV£i)oTg 

Rj  =  ((£j  +  £  J/CjUjOTj 

where  £=£,.+  £  -  £,£  (from  Deardorff  1978) 

1  f  g  f  g 

Hf  =  "PafCpk2z2  S  II  (Balick  et  al.  1981;  Oke  1978) 

Ef  S  (pafCp/0,66)le8(V  *  e(Ta^/(ra  +  rc)  (Gates  198°5  Lee  1978)* 


Ground 


F  =  O  -  af)a  S  ♦  R  -  r  -  «  -  E  -  G 
8  *  8  8*  8*  g  g 


(1  -  o.)ot  S:  o,  and  a  are  inputs;  S  as  above, 
t  8  l  g 


V  S  (l  “  °f)Rs*  *  °fl*goTg  ♦  0  *  eg)efoTfJ’/£i  (Deardorff  1978) 

V  3  (1  "  *  (1  ‘  £g)RsiI  ♦  OfUgOTg  ♦  (1  - 

(Deardorff  1978) 

Mg  3  pagCpCHgU»f(T8  “  Taf*  1978) 

Eg  “  p*gV»fL(,»  ‘  <*af)  (Deardorff  1978) 


*  References  cited  in  this  Appendix  are  more  fully  identified  in  the 
References  section  at  the  end  of  the  main  text. 


Al 


G  =  -K(T  -  T  )/Z  T  is  the  temperature  at  the  first  grid 

X  g  X  X 

point  below  the  surface,  distance  Z  ,  in  TSTM  heat  transfer 
algorithm. 


Additional  Definitions  and  Relationships  Needed 
to  Evaluate  Heat  Budget  Terms 


R  ,  =  oT4(0.61  +  0.050  Ve  )  O  +  K.n2) 
s+  a  3  i 

where  n  is  the  decimal  portion  of  cloud  cover  and  is  a 

coefficient  dependent  on  cloud  type  (Balick  et  al.  1981;  Oke 
1978;  Sellers  1965). 

AZ  =  Z;  Z  is  the  height  above  the  ground  where  the  temperature 
and  wind  speed  are  measured. 

AU  =  U  -  U  r;  U,  is  defined  below, 
a  af’  af 

A8  a  6(Ta)  -  0(Tf)  =  6(Ta  -  Tf)  where 

6(T)  =  T(1O00/P)0*2®^  ;  definition  of  potential  temperature. 

p  t  =  0. 3481*/ 1 (T  ♦  Tf)/2)  ;  from  the  ideal  gas  law  assuming  dry  air 
(Lee  1978)* 

es(Tf)  =  6.108  exp)17.269(Tf  -  273.16)/(Tf  -  35.8))  Teten's  equation 
for  saturation  vapor  pressure  as  given  by  Hurray  (1967). 

e(T  )  =  e  (T  )/RH  where  Ml  is  the  relative  humidity, 
a  s  a 

2 

r  -  U  /lit  is  the  friction  velocity  defined  by 

U*  a  kUa/Un|(2a  -  Zd)/Z0)r)  (Oke  1978) 

0  979 

where  2^  a  0.701  *  2^'  is  the  sere  displacement  height 
and 

o  997 

2fl  a  0.131  *  2^*  is  the  roughness  length.  2^  is  the 
foliage  height  (Rosenberg  1974). 

I'  a  Cj  (1  -  C^Hi)  is  a  correction  factor  for  stability 
(Oke  1978). 


M 


Cj  =  l.o  C2  =  5.0  C3  =  2.0  if  Ri  >  0  (Oke  1978) 

C  =  1.175  C,  =  15.0  C3  =  0.75  if  Ri  <  0 
(Lamb  1974) 

Ri  =  {g/[6(Ta;  +  6CTf) ]/2} [(AG/AZ)/(AU/AZ)2] 

where  g  is  the  acceleration  due  to  gravity.  If 
Ri  >  0.2,  Ri  is  set  equal  to  0.199. 

ra  =  {ln[  (Z  -  Zd)/Zo]H2/(k2Ua). 

rc  =  /1AI  and  is  the  canopy  resistance  to  water  vapor  diffusion 
where 

r  =  (0.05  +  0.0021S)  1  (Gates  1980)  and  is  the  stomatal 
s 

resistance 

LAI  -  7o^  (Deardorff  1978)  and  is  the  leaf  area  index. 

p  =  0.348P/T  .  as  for  p  - 
rag  af  af 

T  _  =  (1  -  ar)T  +  ct(0.3T  +  0.6Tr  +  0.1T  )  a  mixture  equation 
af  f  a  a  f  g 

(Deardorff  1978). 

Ctr  =  (1  -  ar)Cu  +  arC„  ;  heat  transfer  coefficient  equation, 

Hg  f  H  f  H. 

°  oh 

interpolated  between  ground  with  no  cover  (C^  )  and  complete 

o 

cover  (Cu  ) (Deardorff  1978) 

Hh 

Cu  =  k2/ln(Z  /Z  )2 
H  '  a'  u 
o 

C  =  k2/iln[(z  -  Zd)/Z0l)2  . 
h 

U  £  =  0.83O.U  cV2  +  (1  -  oJU  ;  mixture  (Deardorff  1978). 
at  t  a  l  a 

<lf  =  r"qs(Tf)  +  (1  -  r")qflf  (Deardorff  1978) 
r"  =  ra/(r8  +  r.) 
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qs(Tf)  =  0.622/lP/e  (Tf)  -  0.378]  . 
q(T  )  =  0.622/[P/e(T  )  -  0.378]  . 

qaf  =  (1  ‘  af)q(Ta)  +  aff°-3^(Ta)  +  0,6qf  +  0>lqgJ 
qg  =  Wqs(Tg)  +  (1  -  W)qaf  . 

Basic  Symbols 

Energy  budget  component  -  flux  per  unit  area 

S  Solar  irradiance  at  the  top  of  the  canopy 
R  Thermal  infrared  irradiance 
H  Sensible  heat  exchange  with  the  atmosphere 
E  Latent  heat  loss  to  the  atmosphere 
G  Conduction  of  heat  in  the  top  soil  layer 

Physical  properties 
T  Temperature 
P  Pressure 
p  Air  density 
U  Wind  speed 
q  Specific  humidity 

Coefficients 

of  Effective  foliage  cover,  decimal  fraction 
c  Graybody  thermal  infrared  emissivity 
or  Solar  absorptivity  (1  -  reflectivity) 

C  Heat  transfer  coefficient 

W  Relative  saturation  of  the  ground  near  the  surface  with 

respect  to  field  capacity  (0  -  1.0) 

X  An  arbitrary  multiplier  (x  >  0)  of  r  used  to  account  for 

s 

such  factors  as  senescence,  stress,  etc. 
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Basic  subscripts 

a,f,g  Air  (at  instrument  height),  foliage,  ground  surface 
af  Air  within  foliage  layer 
ag  Air  near  the  ground  surface 
I , t  Downwelling,  upwelling 

n  Net  exchange  (thermal  infrared)  between  foliage  and  the 
ground 

s  Sky  (Rg)  or  saturation  (es»<ls) 

c  Canopy  (ensemble  of  foliage  elements)  value  for  resistance 
to  water  vapor  difTusion 

Physical  properties 

Specific  heat  at  constant  pressure  of  dry  air 
k  von  Karmon  constant  (0.4  used  here) 

K  Heat  conductivity  of  the  soil  (assumed  constant) 

L  Latent  heat  of  evaporation  L  =  f(T  ) 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Balick,  Lee  K._ 

Inclusion  of  a  simple  vegetation  layer  in  terrain 
temperature  models  for  thermal  infrared  (IR)  signature 
prediction  :  final  report  /  by  Lee  K.  Balick,  Randy  K. 
Scoggins,  Lewis  E.  Link,  Jr.  (Environmental  Laboratory, 
U.S.  Army  Engineer  Waterways  Experiment  Station).  — 
Vicksburg,  Miss.  :  The  Station  ;  Springfield,  Va.  : 
available  from  NTIS,  [1981], 

3^,  (53  p.  :  ill.  ;  27  cm.  —  (Miscellaneous 
paper  /  U.S.  Army  Engineer  Waterways  Experiment 
Station  j  EL-81-1*) 
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tinder  Project  No.  ltA762730ATlt2 ,  Task  Ah,  Work 
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1.  Computer  simulation.  2.  Relief  models. 

3.  Infra-red  detectors,  it.  Temperature.  5.  Vegetation 
mapping.  I.  Scoggins,  Randy  K.  II.  Link,  Lewis  E. , 
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